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A distinctive feature of layered conductors is the presence of low-energy electronic collective 
modes of the conduction electrons. This affects the dynamic screening properties of the Coulomb 
interaction in a layered material. We study the consequences of the existence of these collective 
modes for superconductivity. General equations for the superconducting order parameter are de- 
rived within the strong-coupling phonon-plasmon scheme that account for the screened Coulomb 
interaction. Specifically, we calculate the superconducting critical temperature T c taking into ac- 
count the full temperature, frequency and wave-vector dependence of the dielectric function. We 
show that low-energy plasmons may contribute constructively to superconductivity. Three classes of 
layered superconductors are discussed within our model: metal-intercalated halide nitrides, layered 
organic materials and high-T c oxides. In particular, we demonstrate that the plasmon contribu- 
tion (electronic mechanism) is dominant in the first class of layered materials. The theory shows 
that the description of so-called "quasi-two-dimensional superconductors" cannot be reduced to a 
purely 2D model, as commonly assumed. While the transport properties are strongly anisotropic, it 
remains essential to take into account the screened interl&yer Coulomb interaction to describe the 
superconducting state of layered materials. 



PACS numbers: 

I. INTRODUCTION 

Recent years have witnessed the discovery of many new 
superconducting materials: high-temperature cuprates, 
fullerides, borocarbides, ruthenates, MgB2, metal- 
intercalated halide nitrides, intercalated Na^CoOai, etc. 
Systems such as organics, heavy fcrmions, nanoparticles 
have also been intensively studied. Many of these systems 
belong to the family of layered conductors characterized, 
e.g., by strongly anisotropic electronic transport proper- 
ties. Recently, it was reported that even Bai-^K^BiOs 
has a layered structured An interesting question raised 
by the observation of superconductivity in all the sys- 
tems mentioned above is the following: why is layering a 
favorable factor for superconductivity? The present pa- 
per addresses this question. We show that layering leads 
to peculiar dynamic screening of the Coulomb interac- 
tion and that this is important for the description of the 
superconducting state in layered conductors. 

The conventional theory of superconductivity has 
mostly dealt with three-dimensional (3D) isotropic sys- 
tems, although some papers have also described the im- 
pact of the Fermi surface anisotropy on the supercon- 
ducting state (see, e.g., review Ref.yJ). In this theory the 
Coulomb repulsion is described by a static pseudopoten- 
tial \i* and its value is reduced because of the well-known 
logarithmic factor \n(E/H.) where E is an electronic en- 
ergy and f2 is a characteristic bosonic (e.g. phonon) en- 
ergy. Such a static approach is justified by the large value 



of the plasmon frequency f2 p z(q = 0) = min{0(q)} = £l p i 
in usual metals, where f2 p z ranges between 5eV and 30eV. 
Such high energies imply a perfect, instantaneous screen- 
ing of the Coulomb interaction. 

Layered conductors have a structure of the plasmon 
spectrum that differs fundamentally from 3D metals. In 
addition to the high energy "optical" collective mode 
mentioned above, the spectrum contains also an impor- 
tant low- frequency part (see below). The screening of 
the Coulomb interaction is incomplete and the dynamic 
nature of the interaction becomes important. As a re- 
sult, the interplay between the attractive interaction and 
the Coulomb term is more subtle than introduced in 
the conventional theory of superconductivity. It is on 
this screened Coulomb term and its interplay with the 
electron-phonon mechanism that we focus in the present 
paper. 

Our goal is to evaluate the additional impact of dy- 
namic screening on pairing in layered superconductors. 
The pure plasmon mechanism (that is, in absence of any 
other attractive interaction) has been discussed previ- 
ously for 3D and 2D systems (see, e.g., Refs. l4l5l lal7l) The 
acoustic plasmons for spatially separated layers in metal- 
oxide semiconductor structures were introduced and an- 
alyzed by Y. Takada in Ref. 0. The author indicated the 
possibility of acoustic-plasmon mediated superconductiv- 
ity. In the present paper we focus on layered conductors. 
More importantly, we consider plasmons' contribution in 
conjunction with the phonon mechanism. It is assumed 
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that the phonons themselves provide the pairing so that 
at T = OK the compound is in the superconducting state. 
In other words, the presence of phonons is sufficient 
to overcome the static Coulomb repulsive interaction. 
Within this scenario the dynamic screening acts as an ad- 
ditional factor. Therefore, in the absence of the plasmon 
term we obtain the conventional Eliashberg equations; 
the electron-phonon coupling constant and the Coulomb 
pscudopotcntial are thus considered as parameters to be 
determined from experimental data (see, e.g., Ref. 0). 
Note that the contribution of phonons and plasmons to 
the superconducting state has aslo been considered in 
Ref. Hol for fullerides. We also point out that we consider 
the electron-phonon interaction (phonon-plasmon mech- 
anism) for concreteness. However, since our attention 
is set on the Coulomb contribution to the total pairing, 
our approach is valid for other mechanisms as well. The 
advantage of the present approach is that we are not re- 
stricted to answer the question whether or not plasmons 
themselves can provide superconductivity, but allows to 
answer the question whether low-energy plasmons can 
sustain or enhance the pairing induced by other mecha- 
nisms. 

We discussed briefly our approach in Refs. ITllH^ITMl^ . 

The present article contains a detailed analysis of the di- 
electric function, the plasmon spectrum and its impact on 
superconductivity in layered superconductors. Further- 
more, we apply the theory to characteristic examples of 
three classes of materials: the metal-intercalated halidc 
nitrides, organic and high-tcmpcrature superconductors. 

The structure of the paper is as follows. In scctionlTTlwe 
present the main equations describing layered supercon- 
ductors and discuss the electron-phonon and Coulomb 
contributions to the pairing-interaction kernel. In section 
rTllwe discuss the dynamic screening of the Coulomb in- 
teraction in layered conductors. The dielectric function 
and the resulting electronic collective excitations (layer 
plasmons) will be described. It is essential that the di- 
electric function is evaluated and analyzed in the ther- 
modynamic Green's function formalism; this allows us 
to calculate T c . In the next section, Sec. IIVI we con- 
sider three classes of layered superconductors: the metal- 
intercalated halide nitrides, the organics and high-T c su- 
perconductors. The conclusions are presented in Sec. |VJ 



II. MAIN EQUATIONS 

We consider a layered system consisting of stacks of 
conducting sheets along the z-axis separated by dielec- 



tric spacers (see Fig. Because of the conductivity's 
high anisotropy it is a good approximation to neglect 
transport between the layers (see Sec. lIV|l . On the other 
hand, the Coulomb interaction between charge-carriers is 
effective both within and between the sheets. To ensure 
charge-neutrality we further introduce positive counter- 
charges spread out homogeneously over the sheets. 
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FIG. 1: The layered electron gas (LEG) model. The con- 
ducting sheets (dark) are stacked along c and separated by 
spacers (light) with dielectric constant tu- The model con- 
siders an infinite stacking of layers. The electrons are moving 
within the conducting sheets. The Coulomb interaction is ef- 
fective within, but also between the sheets (see text). L is the 
interlayer distance. 



In order to calculate the critical temperature T c for the 
superconducting transition we start with the equation for 
the superconducting order parameter A(k, u> n ) and the 
renormalization function Z(k,LU n ): 
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T 



(la) 

d 3 k' 

^— ^r(k, k';o>„ - u} m )F f (k',uj m ), 

(lb) 

22 J ^p-r(k,k';w„ - w m )G(k',w m ), 



where =< |C^ k ^ > is the Gor'kov pairing function, 

G =< g.Ck.o- > is the usual Green function, and T is 
the total interaction kernel. 

These equations can be rewritten in the following form 
(at T = T c ) 
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In these equations and in the rest of the paper we use 
the following notations: k = (k||,fc z ), where the z-axis 
is chosen to be perpendicular to the layers. We use 
the thermodynamic Green's function formalism (see, e.g., 
Ref. Il5|) with u„ = (2n + l)vrT. Because of the re- 
lation u n — uj m = 2itT{n — m) wc often use the short 
hand (n — m) to denote the frequency dependence [e.g., 
(n + m + 1) stands for to n — u>_r m +n]. Finally, we define 
n (k) = A„(k)Z„(k), w„(k) ee w n Z n (k), A n = A(w n ) 
and Z n = Z{uj n ). 

For concreteness we focus on the case where the inter- 
action kernel is a sum of electron-phonon and Coulomb 
interactions. Then, the total kernel T = r(q, uu n — w TO ), 
with q = k — k', is written in the form 
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D(q,n — m) is the phonon Green function, fi^(q) the 
phonon dispersion relation; summation over phonon 
branches v is assumed. The second, important Coulomb 
term r c is written in its most general form as the ratio 
of the bare Coulomb interaction V c (q) and the dielectric 
function e(q,w n — w m ). Both the Coulomb interaction 
and the dielectric function have to be calculated for a 
layered structure. It should be noted that the relation 
between these two interactions [r p ^ and r c , Eq. (J2J] is 
more subtle. For example, the Coulomb screening af- 
fects the value of the electron-phonon matrix elements 
(see, e.g., Ref . ITfit) . Here, however, we do not calculate 
the electron-phonon coupling constant A [see Eq. (|19a|l ] 
and, similarly to the treatment of conventional supercon- 
ductors (see, e.g., Ref. 0), we use the values determined 
from experimental data. For example, in the case halidc 
nitrides considered in Sec. IIVBI A was determined from 
heat capacity measurements^. 

The Coulomb potential V c (q) is the Fourier transform 
of the 3D Coulomb interaction V c (r) = e 2 /ej\/|r|, where 
6m is the dielectric constant of the spacers, and takes the 
form (see appendix^) 



Vc(q) 
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R{q\\,Qz), 



(6) 



where R(q\\,q z ) is defined in Eq. JHJ below. Introduc- 
ing dimensionless quantities q = q\\/2kp, kf = 2kpL (L 
is the interlayer distance, kp the in-plane Fermi wave- 
vector) as well as N(0) = mt,/2Trh , the 2D electronic 
density of states (mj, is the band mass), we can write 



K(q) = 



^ C R(q;g z ) 

N(0) q ' 



with 



R{q, qz) 



sinh(/Ci?g) 
cosh.{npq) — cos(q z L) 



(7) 



(1 - <*q.0). (8) 



Eq. JO} contains the product of the two-dimensional 
Coulomb potential and the function R(q\\,q z ) which re- 
flects the layered nature of the studied system. As ex- 
pected, limi^oo i?(q) = 1, whereas i?(q) = 2<7||/|q| 2 L for 



L< 1. Furthermore, V c diverges as l/|q| for |q| 



0. 



in agreement with the 3D character of this limit. Note 
that R(q; q z ) contains the factor (1 — 5 qi o) reflecting the 
presence of the neutralizing positive ion counter-charges; 
the presence of this term is implicit in the following. 

Eq. J7J contains the dimensionless Coulomb interaction 
constant defined by 
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vp is the Fermi velocity, c the vacuum speed of light and 
a is the fine structure constant. Note that A c ~ r s , where 
r s = \f2jkprB {tb = h 2 eM/me 2 is the Bohr radius) 
is the well-known dimensionless electron density radius 
defined here for a layered electron gas. 

The electronic screening of the Coulomb interaction is 
described by the dielectric function e(q, uj n ) written in its 
most general form as 



e(q,w„ -u) m ) = 1- V r c (q)n(q,w Ii 



(10) 



og we 

(see, e.g., Ref. Il7j) for real 3D metals the RPA provides 
a qualitative description, whereas a quantitative analy- 
sis requires to go beyond this approximation. For the 
systems of interest RPA is favorable because of the in- 
equality A c < 1 (see below). Note that the contribution 
of the background dielectric function and the inequality 
cm > 1 could be essential (cf. e.g. Ref. ll9T) . It would 
be interesting to perform more exact calculations using 



4 



methods as those of Ref. |2CJ and, in addition, take into 
account the band structure of real materials instead of 
the LEG model. We think that the approximation based 
on the inequality A c < 1 provides the adequate physical 
picture. 



III. LAYERED CONDUCTORS: ELECTRONIC 
COLLECTIVE MODES 

A. Plasmon bands 

The spectrum of collective electronic excitations is de- 
termined by the poles of the two-particle Green function 
which coincides with the poles of the vertex r c 
The latter is the analytic continuation (see, e.g., Ref.| 
of the function r c (q, w„), Eq. ©. These poles corre- 
spond to the zeros of the real-frequency dielectric func- 
tion 1 - y c (q)n(q||, uj) = 0. At T = the real part 
of the polarizability of a single layer takes the form 
(ui > kq^vp) 



(o/n„, 



(q,<£l- 

cf. 15) 
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Re{IL(q,uj)} = 2N(0) 
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^Juj 2 - {hq\\v F y 
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For uj 3> hquVF this expression reduces to Re {IT} ~ 
N(Q)h 2 q 2 v 2 F /uj 2 , as obtained in Ref.H3 

From Eqs. Q, (JSJ) and we derive the general ex- 
pression for the plasmon dispersion relation: 



Kq\\v F > 



(N(0)V c ) 2 
\ + N(0)V c 



(12) 



where V c = V c (q,q z ) is the Coulomb interaction defined 
in Eq. Q. If N(0)V c » 1 we obtain the optical plas- 
mon uj = hquVp-^/l + N(Q)V C (this corresponds to the 
hydrodynamic approximation for small q; see Ref. l2l|) . 
The plasmon band uj = ui(qn,q z ) is confined between 
the upper branch with q z = (in-phase motion of the 
charge carriers) and the lower branch at q z = tt/L (out- 
of-phase motion of carriers). Indeed, for uj 3> fiq\\VF 
Eq. (jT2")l reduces to the expression uj ~ Tiq\\VF y/N(0)V c 
which at q z = leads to the usual "optical" plasmon with 
Vt 2 pl = uj 2 { q]i] = 0,q z = 0) = Ae 2 e F /e M L. For q z = tt/L, 
on the other hand, we obtain the dispersion law for the 
"acoustic" plasmon (linear in q\\ ) of the form of Eq. Q12JI 
with 



N(0)V c (q,q z = T 



A c smh(Kpq) 
q COsh(Kpq) 



1 



X c kplil3) 



For q z = tt/L and q»L <C 1, we obtain uj ks (Cl p iL/2)qu. 

Thus, the plasmon spectrum of a layered conductor, 
Eq. (121) has the rather complicated structure shown in 
Fig.H 




q,/2k F 

FIG. 2: Electronic excitation spectrum for the layered elec- 
tron gas (LEG) (see also Ref. Il3h . The solid, dotted and 
dash-dotted lines in region I (uj > hqnVF) are plasmon dis- 
persion relations corresponding to the q z indicated on the fig- 
ure. The area delimited by the q z — (solid) and q z — tt/L 
(dash-dotted) bands contains the dispersion relations for all 
q z . The branch q z = tt/5L (dotted) is shown as an exam- 
ple. The hashed area denotes the electron-hole excitation 
continuum in which plasmon (Landau) damping occurs. The 
long-dashed line separates region I (uj > hq^vp) from region 
II (uj < 7m7||«f). 



The dispersion can be viewed as a continuous set of 
"acoustic" branches parametrized by q z G [Q,ir/L] (the 
slope of the acoustic plasmon at 511 — > is smallest for 
q z = tt/L and increases as q z — > 0). Only the upper 
branch for q z = represents an "optical" branch and, as 
expected for the long wavelength limit, corresponds to 
the usual 3D plasmon. Crucial for the phenomenon of 
dynamic screening and its impact on the pairing is the 
presence of the low-energy collective excitations, which 
can play a role similar to phonons (they can be labeled 
"electronic" sound). 

Note that the low-energy plasmon branches, so-called 
"demons" , also appear in the presence of different over- 
lapping bands (e.g. "light" and "heavy" carriers; see, 
Ref. and the review Ref. 0). We emphasize that the 
case considered in the present paper is entirely different. 
Indeed, the appearance of "acoustic" branches is caused 
by the presence of spatially separated conducting layers 
and the out-of-phase motion of the carriers in neighbour- 
ing planes. 

The partial density of states can be determined for 
each plasmon band (corresponding to each q z ) from the 
dispersion relation, Eq. I] 12(1. As first pointed out by two 
of the authors, Bozovic et al. in Ref. 0, the density of 
states considered as a function of energy is peaked at the 
boundaries, that is for q z = and q z = tt/L. A good 
approximation is thus to model the plasmon spectrum of 
a layered conductor as consisting of two branches: the 
upper "optical" branch (q z = 0) and the lower "acous- 
tic" branch (q z = tt/L). We have shown earlier— that 
the optical branch gives an essentially repulsive contri- 
bution to the pairing interaction. We therefore include 
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this latter part in the effective repulsive /x* of region II 
(Fig. [2J • In the following we consider only the contribu- 
tion of the dominant acoustic branch at q z = tt/L. The 
contribution of all other branches only enhances the effect 
of the q z = tt/L branch, as we discuss below. It is worth 
emphasizing that the existence of the latter branches is 
specific to layered materials. 

We end this section by noting that the inclusion of a 
residual interlaycr hopping would imply the appearence 
of a small gap for the acoustic plasmons. The size of this 
gap is determined by the intcrplane hopping parameter 
t z . The more isotropic the system becomes, the larger is 
the gap. In the isotropic limit the one degenerate opti- 
cal plasmon branch observed in 3D metals is recovered. 
As mentioned in the introduction, however, the materi- 
als of interest for the present paper fSec. HV|) have a ratio 
t z /tn < 10 -3 so that discarding interlaycr transport is a 
good approximation. Fu rther support for this approx- 
imation is found in Ref. I22l28l (and references therein) 
from dielectric properties and lattice dynamics studies of 
high-temperature superconductors. 

B. Screening of the Coulomb interaction: the 
dielectric function at finite temperature. 

To study the impact of dynamic screening on the 
superconducting state we need to calculate the dielec- 
tric function, Eq. 1|10[) . which contains the polarizabihty 
II(q, u> n ). In particular, to obtain T c , we have to deter- 
mine these functions at finite temperatures. In RPA the 
polarizabihty takes the well-known form 

n (q|l ,^) = 2 / d% . ^ " /k " +q " . (w) 

J «^„ + £ k|| -£k||+q|| 

where fk = /(^ky ) is the Fermi distribution and all wave 
vectors lie in the plane of the layered structure. To the 
best of our knowledge all previous works concerned with 
layered structures were done either us ing, the calculated 
polarizabihty at T = (see, e.g., Ref. |2J) or taking the 
static limit for non-zero temperatures (as done, e.g., in 
Ref Hi Here we calculate the polarizabihty both at fi- 
nite temperatures (using the temperature Green's func- 
tion formalism) and all values of (q,w„). As the cal- 
culations of the next sections will show, the tempera- 
ture dependence of the polarizabihty can be neglected in 
some cases (e.g. for halide nitrides, Sec. IIV B)l but should 
be taken into account for the cuprates where the ration 
T c /ep is not negligibly small (see Sec. IIV D|) . In gen- 
eral, the proper account of dynamic screening requires to 
consider all three parameters. Note that to render the 
numerical problem tractable when solving Eq. (|1 7|) be- 
low, we reduce the number of integrals to be performed 
numerically by writing Eq. (|14|) in polar coordinates and 
integrating analytically over the angles (cf. appendix [5]). 
The remaining /c-integration is then done numerically 

Fig. 01 displays the polarizabihty, Eq. (|14f) . of the 
electron gas of a layer, as a function of wave-vector 



q = q\\/2kp for different values of the frequency u n and a 
typical temperature T/ep = 0.03 which applies to high- 
T c superconductors (see Sec. IIVD(I . We first point out 
that the g-dependence of this function is essentially re- 
stricted to the interval [0, 4fcf]- 
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FIG. 3: Electronic polarizabihty in RPA as a function of 
q — 9||/2/cf for various values of w n — Um = 2(n — m)nT 
(from top to bottom) and a typical temperature T/ef = 0.03. 
The lowest, dash-dotted curve is obtained for T = OK and 
remains the same for all frequencies [i.e. at T = 0, Tl(q < 
2k F ,uj n -uj m ) =-2N(0)]. 



Let us now consider the temperature dependence of 
the polarizabihty. First we discuss the case T = 0. For 
the real-frequency polarizabihty the analytical continu- 
ation gives the result derived by Stern in Ref. OOJ. On 
the other hand, in the Matsubara temperature-Green's 
function formalism the polarizabihty does not depend 
on frequency for uj < quVp (cf. Ref. Il5k . Sec. 20.2): 
II(<7, u> n ) = — 2A(0) = —nib/nTi 2 (bottom dotted line of 
Fig. |3J). For T > the polarizabihty vanishes at high 
frequencies as shown in Fig. 0J The higher the tem- 
perature, the smaller is the frequency-range over which 
the polarizabihty remains finite. Note that the shape of 
II(<7, u> n — uj m = 0) (lower solid line on the figure) is al- 
most unaltered until very high temperatures. This can 
also be directly seen from the analytical expression [see 
appendix, Eq. (|B8|) ]. 



Using the above results for the polarizabihty we calcu- 
late the dynamic dielectric function, Eq. I|10|) . or rather 
it inverse e _1 (q,w n — oj m ) since it is this quantity that 
enters the vertex r c , Eq. JSJ. The result is shown on 
Fig. 0| for the same values of parameters as in Fig. 01 We 
point out a few important properties of the inverse di- 
electric function. This function is bounded for all q and 
uj n — uj m . For high frequencies and/or large wave- vectors 
e -1 (q, ui n — uj m ) — > 1 , meaning that the Coulomb interac- 
tion is unscreened in these cases. In fact, it can be shown 
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FIG. 4: Inverse dielectric function of the LEG. The solid, 
dashed, and dotted lines correspond to q z L = 0, n/2, n respec- 
tively. The three bottom curves are obtained for u) n — ui m = 0, 
the three middle curves are for n — m = 10 and the three up- 
permost curves are for n—m = 40. The other parameters were 
given in the previous figure. Note that e~ 1 (q — > 0,cu n — ^m) 
is zero (perfect screening) only in the static limit. 



that for Tiqvp <C Lu n 
the formii 



the dielectric function takes 



e(hqv F <C uj n 



j) ~ e D (n ~ m), 



where 



coin — m) = 1 + 



n 2 



(15) 



(16) 



Eq. (|15|) describes the dielectric function in the Drude 
limit. Note that latter expression is exact in the limit 
q = and lu > 0i£ 

The result of Fig. 0] shows that for any finite fre- 
quency, the long wave-length limit takes the form e(|q| — > 
0,LJ n - uj m > 0) — > £_d(w„ - w TO ). Thus, it is only in 
the static case uj n — u> m = that the Coulomb inter- 
action is "perfectly" screened (exponential screening in 
real space): e(q — > 0, u) n — io m = 0) « [1 + Ky F /|q| 2 ] -1 . 
The latter limit is the so-called Thomas-Fermi screening 
of the Coulomb potential. For the LEG the screening 
length is given by n^p = nN^Q^. In all other cases 
the limit of long wave-lengths is given by the Drude limit. 

The dielectric function describing the screening in lay- 
ered conductors will be used in the next section to cal- 
culate the effect of the dynamically screened Coulomb 
interaction on T c in several classes of layered supercon- 
ductors. We thereby use the full wave-vector, frequency 
and temperature dependence of e(q, uj n — uj m ) calculated 
in this section. 



IV. APPLICATION TO VARIOUS LAYERED 
SYSTEMS 

In this section we consider the phonon-plasmon mecha- 
nism and in particular the impact of the dynamic screen- 
ing of the Coulomb interaction on the superconducting 



state of several layered systems. We discuss specific ex- 
amples belonging to three classes of materials: metal- 
lochloronitrides, organics and high-T c superconductors. 
To this aim, we first rewrite Eqs. (|2al2b|) in a form ad- 
equate for layered conductors and convenient for calcu- 
lations. We then evaluate the critical temperature T c of 
the various compounds. 



A. Numerical analysis 

We assume isotropy within the layers. The order pa- 
rameter and the renormalization function can therefore 
be written as A„(k) ~ A n (k|i = kp,k z ) = A n (k z ) and 
Z„(k) = Z n (k z ) (see appcndixQ . This approximation is 
valid since the Cooper instability (see, e.g., Ref. fl5^ and, 
correspondingly, the pairing, occurs on the Fermi surface. 
As for the integrands in Eqs. 1)1150 . they depend mainly 
on the momentum transfer . Thus, in the layered elec- 
tron gas the order parameter depends on the frequency 
and the wave-vector perpendicular to the layers. 

To perform numerical calculations with Eqs. I|2al2b(l 
we follow the standard procedure adapted to the case of 
layered materials. We first express the integral over ky 
in terms of an integral over energy and in-plane wave- 
vector amplitude and carry out the former analytically 
(see appendix Q. We thereby reduce the equations to a 
form containing one-dimensional k-space integrals. Note, 
however, that the dielectric function also contains an in- 
tegral to be performed at each iteration of the calculation 
[Eqs. (|10I14|| ]. The resulting equations take the form (ap- 
pendix Q 
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where N z are the number of k z points taken in the first 
Brillouin zone and the kernel is given by 



r(<7z;"~ m ) = AZ?(n — to) + A c r c (g z , n — to), (18) 



with 



Xsm * M . )|2 _ (19a) 



7T Jo \f\^q" 
T c (q z ,n- m) = (19b) 
1 f 1 dq R{q,q z )/q 



Wo y/l - q 2 e(q,q z ,n- m) 

All other quantities were defined in Sec.^ Two simplifi- 
cations are made in the following calculations that allow 
to single out the effect of low-energy electronic collective 
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modes on superconductivity. The first is to replace the 
phonon contribution to the pairing to one (or two, see be- 
low Sec lIVCjl characteristic phonon modes. The second 
is that we can set k z = n/L, based on the analysis made 
in Sec. lIII Al which shows that this wave- vector gives the 
largest contribution to the pairing. Consequently, the or- 
der parameter and the renormalization function taken at 
the zone boundary along k z are function of frequency u> n 
only. We emphasize, however, that these two simplifica- 
tions are not affecting the main results presented below. 
Rather, they would lead to smaller coupling constants 
necessary to reach a specific critical temperature T c . For 
example, we expect that taking into account all plasmon 
acoustic plasmon branches would lead to higher criti- 
cal temperatures (or correspondingly to smaller Coulomb 
coupling constants for a given T c ) since they contribute 
to the attractive pairing interaction, though with lesser 

weight iiii 

As discussed in Sec. IIII Al the excitation spectrum of 
a layered electron gas (see Fig. |2J allows us to divide 
the (gii , w)-space into two main regions. The first (region 
I) contains the additional collective excitations discussed 
here and this region should be considered exactly in the 
equations above. On the other hand, region II contains 
no such additional plasmon features and the equations 
have a form analog to the 3D case in this area. Therefore, 
we make a further step by dividing the Coulomb part, 
Eq. (|19bJ) . into two parts corresponding to the two regions 
of Fig. |2 

T c (q z ,n - m) = T l c (q Zl n- m) + T 1 J(q Zl n- m) = (20) 

1 f /"*" + f 1 ) dq R(q,q z )/q 
7T [Jo Jq c J v 7 ! - <f e (^' qz,n-m)' 

where q c = min{2fci?, \oj n — uj m \/hvF} and q c = q c /2kp. 
The part T\ will be considered exactly, in particular with 
respect to the frequency dependence that has the distinc- 
tive features of layered conductors. The part r* 1 , on the 
other hand, is reduced to an effective constant Coulomb 
pseudopotential \i*B(uj n — Q c ) with a standard cutoff 
given by f2 c ~ 10 x 51 (il is the characteristic phonon 
energy). This treatment of region II calls for a comment. 
As was mentioned earlier we consider usual phonon- 
mcdiatcd superconductivity (Eliashberg equations) as a 
starting point of our analysis. Accordingly, the electron- 
phonon coupling constant and the static term fx* are 
treated in the conventional way as parameters to be 
determined from experimental data. Thus, the static 
term fx* is here a phenomenological parameter. We focus 
on the term describing the contribution of the dynamic 
screening. This part (present in region I) will be eval- 
uated explicitely for different systems with the use of 
normal-state parameters (as, e.g., vf or 6m; sec below). 
Note that for the pure plasmon mechanism both static 
and dynamic terms were calculated in Ref. (in 3D). 
This step was crucial since its value was directly related 
to the criterion for the appearence of superconductivity. 
For the phonon-plasmon mechanism, on the other hand, 



we assume that the phonons are sufficient for the oc- 
curence of the superconducting state, which allows to use 
the conventional approach for the static term. Naturally, 
it would be of great interest to calculate the static term. 
Such a full self-consistent calculation including also real 
band structures will be carried out elsewhere. 

As shown in appendix [0 (see also Ref. |!3), the above 
equations can be mapped onto an eigenvalue problem 
written in tensorial form 

K $ = r?$, (21) 

where ®„„ 2 = $ n (n z ) = A n {n z )/y/2n + 1 and K is 
given by Eq. (U3J in appendix Eq. is written 
explicitely in Eq. (|C4(1 . Note that an artificial eigenvalue 
r\ has been introduced in Eq. (121(1 . T c is reached when r\ 
is one. Since all eigenvalues satisfy the inequality 77 < 1 
we only need studying the highest of them. Furthermore, 
the solution of these equations give also the renormalized 
order parameter near T c . We can thus analyse how 
this function is affected by the contribution of low-energy 
plasmons to the pairing interaction. A typical function 
is shown in Fig. [3] and will be discussed in the next 
section. Using the eigenvalue equation (|2T)l [or l|C4|l in 
appendix [C] we apply the theory to various layered su- 
perconductors and calculate their T c . 

B. Intercalated metal halide nitrides 

The first class of materials we consider is the family of 
layered metal-intercalated halide nitrides. We give spe- 
cial attention to this family because low-energy plasmons 
not only contribute to the pairing but, in fact, play the 
key role for the superconducting state, as we show below. 
We believe that this is the first observed system where 
the superconducting state of the electrons is essentially 
self-supported, that is, where the pairing is provided by 
collective excitations of the same carriers as those form- 
ing pairs. 

This family of novel superconductors has been discov- 
ered recently and studied in detail in Refs. Il3l26ll27y 2 ^ . 
I2lll3^ll3lll3ll3ll3ll3lll3ri As is known (see, e.g., Refill) 
the intercalation of alkali atoms and organic molecules 
into the parent compound (Zr,Hf)NCl leads to a super- 
conductor with rather high critical temperature (T c ~ 
25K). Based on experimental 18 i 26 i 27 i 28 i 29 i 30 i 31 i 32 i 33 stud- 
ies and band structure calculations 3 ^ it was concluded 
that electron-phonon mediated pairing is insufficient to 
explain the observed T c , since the electron-phonon cou- 
pling constant appears to be too small. Note also that a 
small nitrogen isotope effect of T c has been observed^ 
In addition, the compounds do not contain any magnetic 
ions and no sign of magnetism has been found in band 
structure calculations. This excludes a magnetic pairing 
mechanism. Finally, normal-st ate proper ties of the mate- 
rials studied intensively in Refs. Il8l35l36l can be described 
by Fermi liquid theory, so that there is no evidence for 
the presence of strong correlations. 
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We apply our approach to this novel layered system. 
We note from Eqs. I|9I17I20[) . that the evaluation of T c 
for a specific compound needs the knowledge of following 
parameters: the interlayer distance L, the band mass mj, 
and Fermi velocity vf, and the dielectric constant of the 
spacers cm- In addition, the evaluation of the phonon 
contribution to the pairing requires the knowledge of the 
characteristic phonon frequency il, the electron-phonon 
coupling constant A and the Coulomb pseudopotcntial 

Specifically we consider Li . 48 (THF) y HfNCl 
(THF=tetrahydrofuranc) as an example, since the 
largest amount of information necessary for the deter- 
mination of T c is available for this material, both from 
experim ent and b and structure calculations. According 
to Refs. EH 0l34l the interlayer distance L and the char- 
acteristic phonon frequency f2 are equal to: L = 18.72A 
and f2 = 60meV. The values of the band mass and 
Fermi energy have been evaluated from band structure 
calculations, Ref. l35i Accordingly to/, = 0.6m e , where 
m e is the free electron mass and £f — leV. For em 
we have chosen the reasonable value cm = 1.75. It 
follows that A c ~ 0.8 and, correspondingly, using 
Eq. JSJ, r s ~ 2 (i.e. close to the high-density limit). The 
value of the electron-phonon coupling constant can be 
estimated from the knowledge of the electron specific 
heat constant 7 and the band density of states N(0). 
Indeed, the electron-phonon interaction renormalizes 7 
as 7 = 7b(l + A), where 7;, = 27r 2 A^(0)/3 is the free 
electron Sommerfeld constant. The value of 7 was 
estimated in Rcf. to be 7 ~ 1.1 mJ/molK 2 , whereas 
band structure calculations^ give N(0) ~ 0.74 eV _1 . 
Thus, A ~ 0.25. Setting u* = 0.1 and using Eqs. 
we obtain T c ~ 24. 5K. The calculated T c is very close 
to the observed value T^ xp = 25. 5K. The essential point 
to note is that in absence of the plasmon contribution 
we obtain TP honon < 1K(!). This demonstrates that, 
indeed, the low-energy plasmon contribution plays a key 
role for superconductivity in metallochloronitrides. 

It would be of great interest to carry out specific tun- 
neling (cf. Ref. |3) and optical measurements on this 
material. We expect that tunneling experiments, simi- 
larly to heat capacity data (see above), will provide the 
value A ~ 0.25, and optical measurements will lead to 
e M 1-75. 



As mentioned earlier, by solving Eq. (|17|) we not only 
obtain T c , but we also get the superconducting order pa- 
rameter $„ = A„/V2n + 1. It is interesting to see how 
<f>„ is affected by the additional pairing arising from the 
presence of acoustic plasmons (Fig. |3J). Note that the 
qualitative form of the order parameter is the same for 
all classes of materials discussed in the present paper. In 
the absence of the plasmon contribution, the order pa- 
rameter is a rapidly decreasing monotonic function of the 
Matsubara frequency (dotted line of Fig. |3J). The effect 
of plasmons reveals itself as an additional "step" in <£> n at 
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FIG. 5: Normalized order parameter $„ = A„/ \J2n + 1 as a 
function of Matsubara frequency u) n = (2n + l)nT. The solid 
and dashed lines obtained for A = 1 and 1.5 respectively, 
represent the order parameter in the presence of the plasmon 
contribution. The dotted line was obtained in absence of the 
plasmon contribution (A c = 0, A = 1). Note the presence of 
an additional structure (step) when including the pairing due 
to low-energy plasmons (solid and dashed lines), fi* — 0.1. 

intermediate frequencies, as examplified by the solid and 
dashed lines on Fig. |SJ It is this positive part of the order 
parameter due to the pairing induced by the low-energy 
collective modes that is responsible for the enhancement 
of the value of T c . We observe that the order parameter 
remains positive over a frequency range also determined 
by the value of A and /i*. Whereas the frequency range 
over which $„ remains positive shrinks with increasing 
A (compare solid and dashed lines on Fig. |3J) it increases 
with increasing fx* (not shown on the figure). This ap- 
parently counter-intuitive behaviour is easily understood 
by the fact that an order parameter extending to higher 
frequencies will pick up more and more repulsive compo- 
nents of the pairing interaction. The "shorter" the step in 
frequency of the plasmon-induced structure, the smaller 
is the repulsive contribution of the effective interaction 
kernel and, consequently, the higher is T c . 

Concluding this section, we emphasize that the dy- 
namic screening of the Coulomb interaction (the con- 
tribution from low-energy electronic collective modes) is 
essential for the understanding of the superconducting 
state in intercalated layered metal halide nitrides. 



C. Layered organic superconductors 

Organic superconductors were predicted in Ref. Wi\ and 
discovered in Ref. |3^. In this section we apply the theory 
to the class of layered organic superconductors (see, e.g., 
Ref. I39ll40ll4lll42ll4^ and references therein) . As an exam- 
ple, we focus on k-(ET) 2 Cu(NCS) 2 (ET=BEDT-TTF is 
a short notation for bisethylencdithiotetrathiofuvalene). 
The basic structural building blocks of these materials are 
large, elongated ET planar molecules stretching along the 
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c-axis forming the conducting layeriaa These thick con- 
ducting layers are separated by thin insulating spacers 
of planar NCS-molecules extending in the <x6-planc. The 
NCS counter-ions take one charge per two ET-molecules 
leaving the ET HOMO partially unfilled. The HOMO 
are 7r-holes delocalized over the large organic molecule 
and form the hole conduction bands. 

Layered organic conductors have highly anisotropic 
transport properties. Typically, the ratio of in- to out- 
of-planc conductivity is at least of the order (t\\/<j± ~ 
-LQ4 | 42 i 4£ 3 an( j structure calculations^ confirm the pres- 
ence of quasi-two-dimensional bands. We emphasize once 
more, however, that only electronic transport proper- 
ties are quasi-2D. As discussed in the previous sections, 
the Coulomb interaction is important in all three dimen- 
sions. In particular, incomplete screening between layers 
implies that carriers from differents layers interact with 
each other, leading to the low-energy electronic collec- 
tive modes discussed here. As we show in the following, 
this aspect is important for understanding the relatively 
high value of the critical temperatures observed in these 
materials. 

Superconductivity has been observed for temperatures 
T < T c ~ 10.4K. Recent studies have shown the im- 
portance of electron-phonon interaction for the pair- 
ing mechanism42i22i5fli5ii52i& For example, isotope effect 
studies of the superconducting T c by isotope substitution 
of C and S atoms on the ET-molecules have singled out 
the effect of intramolecular vibrations for the supercon- 
ducting pairing^2ii& A shift of phonon frequency caused 
by the superconducting transition has also been observed 
with inelastic neutron scattering^ This shift indicates 
that the coupling to miermolecular acoustic phonons con- 
tributes to superconductivity. Further work supporting 
the importance of electron-pho non inter action for super- 
conductivity are given in Refs. I5ll52l53l Therefore, it is 
interesting to apply our phonon-plasmon model to this 
class of materials and study the effect of acoustic plas- 
mons on the superconducting T c . 

As it appears that both inter- and intramolecular vi- 
brational modes are of importance to superconductivity 
we modelize the phonon kernel in Eq. @ by a two-peak 
function 

r p h(u„ - dJ m ) = (22) 

The lower frequency mode £l\ = 5meV corresponds to 
libration and intermolecular modes i^Si The higher fre- 
quency peak is located at the frequency f^, and was 
calculated for the ET- intramolecular vibrations: = 
10meVi£2i£i Given the number of modes present near 
each peak and their possible coupling to the electrons, 
we set w\ = 0.75 and W2 = 0.25. The coupling con- 
stants to each set of modes is then defined as Xj = Xwj 

(i = 1,2). 

To calculate the value of T c , we need to know the value 
of the band mass mj, the interlaycr distance L, the Fermi 



energy ep and the dielectric constant of the spacers cm 
(normal state parameters), as well as the value of the 
electron-phonon coupling constant A and the Coulomb 
pseudopotential fi*. From band structure calculations 
we have mj, = 1.72m e ^ The structure determination 
gives L = 16.2AiiThe average value of the Fermi wave- 
vector obtained from Shubnikov-de Haas measurements 
is kp — 2.6 x 10 7 cm -1 (Ref. 45). Inserting these values 
in ep = h k F /2mb we obtain ep ~ 0.17eV. Note that 
this is exactly the value obtained from band structure 
calculations^ Finally, we extract the value of cm from 
optical reflectance measurements^ 7 - Using Eqs. (|D1ID2|) 
and the data of Ugawa et alwZL we obtain cm = 6.5 (ap- 
pendix [DJ). Note that the ionic screening of the Coulomb 
interaction is more efficient in organics than in metal- 
lochloronitridcs (previous section). One reason for this 
difference is given by the fact that in organics the thick 
conducting slabs are made of large molecules, whereas in 
the metallochloronitrides conducting sheets are thin and 
made of covalently bond atoms. The polarizability of the 
molecules implies better ionic screening of the Coulomb 
interaction and, therefore, a larger value of cm- These 
parameters lead to A c ~ 0.9 and thus r s ~ 2.5. 

The exact value of the electron-phonon coupling con- 
stant A is unknown at present. Estimates range A be- 
tween 0.5 and 1.5m2£2£l Consequently, we present re- 
sults for this range of values in Fig. [S] (we have chosen 
/x* = 0.1 and the cutoff at fi c = 10 x fl 2 = O.leV). Using 
these parameters we calculate T c from Eqs. i|17ll8fl (see 
Fig.©. 
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electron-phonon coupling X 



FIG. 6: T C (A) where A is defined in Eq. 122H . T c is normalized 
to the lowest phonon energy Oi. The upper (lower) curve is 
obtained in the presence (absence) of the acoustic plasmon 
contribution, (i* — 0.1. 



The result shows that the increase of T c in the presence 
of low-energy electronic collective modes is substantial. 
We can quantify this enhancement of T c for the specific 
example studied, the n— (ET)2Cu(NCS)2 compound. Ac- 
cording to our calculation (see Fig. the experimentally 
observed value T c = 10.4 K is obtained for A ~ 1, im- 
plying a coupling to the low- and high-energy phonon 
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modes of Ai = 0.83 and A 2 = 0.28, respectively. Thus, 
k— (ET) 2 Cu(NCS) 2 is an intermediate coupling super- 
conductor. In the absence of acoustic plasmons we obtain 

rpphonon = g 3 K for thig X Thus ^ in pregent cage 40% 

of the value of T c is due to the pairing of electrons via 
the exchange of acoustic plasmons. These calculations 
lead us to conclude that the contribution of low-energy 
electronic collective modes to the pairing is significant in 
organic superconductors (though not dominant as in the 
case of metallohalidenitrides, Sec. II V fj|) . 

D. High-temperature oxides 

In this section we discuss superconductivity in the 
cuprates within our phonon-plasmon model. We ana- 
lyze one specific material, Lai.gsSro.isCuO^, for which 
most parameters have been determined. The normal 
state parameters are: the interlayer distance L = 6.5A, 
the Fermi wave- vector kp = 3.5 x 10 7 cm -1 , the char- 
acteristic phonon frequency Q ~ 15 meV, the dielec- 
tric constant cm — 5 — lOiSiS Therefore, A c ~ 1 and 
r s ~ 2.8. The effective mass m* and electron-phonon 
coupling constant A were determined by Wolf and one of 
the authors from heat capacity measurements The 
obtained values arc A = 2 and m* ps 5m{,. From the re- 
lation m* = (1 + \)m,b we then obtain m& ~ 1.7. Finally, 
the Coulomb pseudopotential is taken to be fi* ~ 0.1. 

The solution of Eqs. (|i7all7bll8) with use of the afore- 
mentioned parameters lead to T c = 36. 5K which is close 
to the experimental value Tf! xp ~ 38 K. It is essential to 
note that in the absence of the screened Coulomb inter- 
action we would obtain T ph = 30 K. Thus, about 20% of 
the observed value of T c is due to "acoustic" plasmons. 
For thin films the stiffness of the lattice usually increases, 
leading to a higher value of the characteristic phonon fre- 
quency f2. Assuming fl = 20 meV, we obtain T c = 49K 
which is close to the experimental value T£ xp = 45 K 
observed, e.g. in Ref. |57|. This value is indeed higher 
than the one of bulk samples. Interestingly, we obtain 
rpph _ 40 g0 that the increase of T c induced by "acous- 
tic" plasmons is again of order of 20%. 

Thus, the dynamically screened interlayer Coulomb 
interaction is important for superconductivity in the 
cuprates. Note that a proper account of the Coulomb 
interaction screening is not only of importance for su- 
perconductivity in these materials, but also for a proper 
description of normal state properties such as lattice 
dynamics 

V. DISCUSSION AND CONCLUSIONS 

The purpose of the present article was to study the im- 
pact of layering on superoenductivity. Particular empha- 
sis was set on the dynamically screened Coulomb interac- 
tion. Layered materials have distinctive low-energy elec- 
tronic collective excitations that provide exchange bosons 



for the pairing between electrons. We have shown that 
these acoustic plasmons lead to an enhancement of the 
superconducting T c . We have applied the theory to three 
classes of layered superconductors: alkali-intercalated 
halide nitrides, organic and high-temperature supercon- 
ductors. 

Within our phonon-plasmon model we observe an in- 
creasing influence of the electronic pairing mechanism 
for the three classes of layered superconductors consid- 
ered. In metal-intercalated halide nitrides the contribu- 
tion arising from low-energy electronic collective modes 
is dominant. These materials arc thus unique in the sense 
that an electronic pairing mechanism is at the origin of 
superconductivity: the exchange bosons are made of the 
same particles (the electrons) than those who bind into 
pairs below T c . In the case of organic layered materi- 
als, the electronic and phononic energies, as well as the 
structure of the conducting layers and insulating spacers, 
leads to a situation where the contribution of phonons 
and acoustic plasmons is of the same order. Finally, in 
the case of high-temperature superconductors, the con- 
tribution of low-energy plasmons is significant but not 
dominant. Within our model the phonon contribution is 
still largest. 

There are other classes of layered superconductors that 
have not been considered in the present article. Among 
them, the most prominent is that of dichalcogcnidcs. We 
believe that some experimental observations^ are related 
to the phenomenon discussed in this article. However, 
we also note that many of the systems belonging to this 
class of materials exhibit charge density wave instabili- 
ties. This both obscures and changes the contribution of 
acoustic plasmons to superconductivity and will be dis- 
cussed elsewhere. 

Another interesting system is the Co02-based layered 
compound studied recently in Ref. We point out that 
the system becomes superconducting only for relatively 
large interlayer distance. This is consistent with the 
present theory and the material deserves further study. 

An essential conclusion of the present work is that the 
physics of layered (super)conductors cannot be reduced 
to the study of one conducting layer (or the layers belong- 
ing to one unit cell as in some high-temperature super- 
conductors). Such simplification relies on the observation 
of "quasi-two-dimensional" transport. However, it misses 
to account for the screening properties of the electron- 
electron Coulomb interaction (and of the electron-phonon 
interaction as well; see Refs. l2^l5^ . As we discussed in 
the paper, the screening is very different in layered ma- 
terials as compared to 2D and 3D isotropic metals. We 
believe that the particular screening properties are essen- 
tial for the behaviour of layered (super)conductors. How 
large the effect of screened Coulomb interaction is, de- 
pends very much on the specific features of the material. 
For example, the covalency within the conducting lay- 
ers and the structure of the spacers, the presence of van 
der Waals gaps will determine its contribution both to 
normal and superconducting state properties. Therefore, 
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the study of screening properties in layered conductors is 
a promising direction to better understand the similari- 
ties and differences between different classes of materials 
and serve as a bridge in the study of properties of 2D and 
3D systems. 
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APPENDIX B: POLARIZ ABILITY 

The RPA polarizability of a single conducting sheet, 
Eq. i|14|) . is written in polar coordinates k|| — (k,ip) 

2 roc 

n(q||,iw„) = -— - / dkkfkI<p{k,q.\\,w n )HBl) 
V Z7T ) Jo 

(B2) 

, ! , ! 1. 

ibJ n + ^ky + Sk||+q|| i^n + Ck|| — ^ky+qy J 

For T = the integral over k|| can be calculated an- 
alytically, leading to the result first derived in Ref. l24l 
At finite temperature there is no simple analytical form. 
However, to reduce the amount of numerical work in solv- 
ing Eqs. P7all7b|) for T c , we calculate the angle integral 
analytically. Using the transformation z = cxp(itp) we 
integrate I v over z in the complex plane to obtain 




h = 



~^[S(k)+9i-9 2 ], (B3) 



APPENDIX A: COULOMB POTENTIAL FOR 
LAYERED SYSTEMS 



Using cylindrical coordinates r = (riuz) (where z is 
perpendicular to the layers) the Fourier transform of the 
3D Coulomb potential V c (r) = e 2 /ej\/|r| is given in lay- 
ered structures by2i 



w = ^E e ~^ ni / dr||e ~ iq 



e 2 e~ 1 
-iqiir,, e e M 



r 2 + (?iL) 2 



-\q z nL 



2ire 2 
1 



-q\\nL 



(Al) 



where we have taken into account the fact that the 
charges are located in the conducting sheets, and thus 
z = nL where L is the intcrlaycr spacing and n indexes 
the layers. Note that the second line of Eq. I|A1[) shows 
how the Coulomb interaction is exponentially decaying 
(in real space) along the direction perpendicular to the 
layers, the decay being determined by q\\L. Performing 
the sum in Eq. ijAlfl we obtain Eq. ©. 

Note that the detailed structure of the spacers separat- 
ing the conducting sheets is not considered in the present 
model. We thus have included the screening resulting 
from polarization effect of the spacers via the dielectric 
constant £m- The dielectric function in the denominator 
of Eq. (JHJ thus accounts for the screening induced by the 
charge-carriers of the conduction bands only. 



with u± = \fA±W, W = VA 2 +B 2 , A = 4fc 2 g 2 (d 2 
C 2 2 -l),B = 4fc 2 q 2 CiC 2 , and 



C = Cl +<2 = ~~ +1-TZ- 

2k 4fcg 



(B4) 



As in the main text we normalize all wave-vectors to 
X = X/2kp, X = q,k. 8j (j = 1, 2) are defined in terms 
of Hcavyside functions as 9j = 9(1 — \z 3 \ 2 ), Zj = Xj + iyj 
with 

Zi = \u+\-Ci, yi = \u-\-(2, (B5) 
%2 = ~\u+\ -Ci, V2 = ~\u-\ - Ca, (B6) 

for B > whereas y\ and yi are interchanged for B < 0. 
Inserting this expression into Eq. i|Bl|) above, we obtain 
the following compact form for the polarizability 

n(« )Wn ) - -jv(o)^ r dAf A ^±[sCk)+9i-m-) 
q 2 Jo 9A 

This expression of the polarizability has been used 
to calculate the dielectric function that appears in 
Eqs. (|17all7b|l and is depicted in Fig. Note that 



lim U(q,u n - ui m ) = -2N(0)f kjF 6 Uni 

q^O 



(B8) 



APPENDIX C: EQUATIONS FOR THE ORDER 
PARAMETER AND THE RENORMALIZATION 
FUNCTION 



We start with Eq. and, as mentioned in Sec. IIV Al 
assume isotropy of the bands within the planes. Thus, 
A and Z depend only on the norm of kii (and on k z ). 
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In this case, it is possible to calculate one of the inte- 
grals over k'n analytically. To this aim, we transform 
the 2D in-plane integration in a way analogous to the 
3D cases generally studied, namely introducing polar co- 
ordinates d 2 k = k'dk' d(p. With fc'dk' = 2ttN(0) d£ and 
q 2 = |k' — k| 2 = k 2 + k' 2 — 2k'kcos(p) the integral over k|| 
is transformed into an integral over energy £ and angle 
p. Using the fact that the resulting energy integrand of 
Eq. (J2J) is falling off as £~ 2 , the main contribution to this 
integral will come from ^/sf *C 1 and we obtain 

d 2 k' = 27rjV(0)d£dy? - 47rN(0)d£ , (CI) 

yi - q 2 

with q < 1. The energy integral can then easily be per- 
formed. Gathering the different terms and assuming that 



A n (k z )Z n (k z ) 



The second equation has been simplified further, reduc- 
ing the sum over m to the range [0, 2n]. The kernel V in 
Eq. l|C2bf) now only contains frequency-dependent terms. 
All frequency-independent terms vanished in the folding 



(q z = k' — k z ) = — = < F( 

U Z ' v/(2n + l)(2m + l) \ 



we finally condense Eq. l|C2Jl to the matrix form (q z = 
k z k z ) 



E Jf z E K ™n{\< - 77 z |) $ m «) = (C4) 
m>0 n' z —l 



n z = l,...,N z . This is the explicit form of Eq. (|2~T|) . 

Note that the kernel K nm (q z ) depends on n and m sep- 
arately and not only on n — m. Furthermore, the kernel 
is even in q z , K nm (q z ) = K nm (\q z \). We have introduced 
the artificial eigenvalue 7/ to map the problem onto an 
eigenvalue equation. T c is obtained when r/ = 1. 



the electron-phonon coupling function (^(q) = g^(q\\), in 
order to define A and r c as in Eq. 1)19(1 , we obtain Eq. 1)17(1 . 
The latter equations have been obtained by discrctising 
k z = -tt/L + 2im z /N z L with n z = 1,...,N Z . Note, 
that the angle-integration can in principle be performed 
exactly, without need of the approximation, Eq. I|C1|1 . 
However, the difference with the present method is min- 
imal and we use the approximation above for simplicity. 



To perform the numerical calculation we cast Eq. or 
(|17|l into a matrix form. We first transform the summa- 
tion over 77i = ... — 1,0,1,... to a sum over non- negative 
7n only. Equation l|17|l then takes the form 



(C2a) 



(C2b) 

I 

of the summation over m. 

Insertin g Eq. (|C2b|) into (|C2a|> . defining < f> n (fc z ) = 
A n (k z )/\/2n + 1 and 



m) + T(q z , n + m + 1) - 8„ tm T'(q z , n - p), > (C3) 

p=o J 

I 

APPENDIX D: DIELECTRIC CONSTANT OF 
THE SPACERS e M 

The dielectric constant of the spacers em can be 
extracted from infrared or reflectivity data. We 
parametrize the dielectric function obtained in these ex- 
periments by the Drude-Lorentz model 

Sus 2 

e{oj) = e 00 +Y J - J 2 3 ■ r + e fc , (Dl) 

J J 

where ef c is the free carrier contribution to the dielectric 
constant. The dielectric constant for the spacers is then 
defined as 

e M = e(u = 0) - e fc = eoo + Sj. (D2) 

j 



- T E W E { f (&, « - m) +f(q z ,n + m +l)}^±, 

rp 2n ^ 7T 

1+7r — Ef E f'(&,n-m). 



m— k'——7 
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For the determination of the dielectric constant of the 
organic material k— (ET) 2 Cu(NCS) 2 we use Ugawa et 
ai.'s reflectivity measurements^ Their parametrization 
gives = 3.2 and 

wi = 0.16 lo 2 = 0.28 uj 3 = 0.47eV, 
O pl = 0.093 n p2 = 0.7 fl p3 = 0.44eV. 



With the correspondance fi^ = SjW? we have Si ~ 
0.762, S 2 ^ 1.581, S 3 ~ 0.968. From these data and 
Eq. JDS} it follows that 

e M = 6.5 (D3) 
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